Multiple Schrodinger pictures and dynamics in shortcuts to adiabaticity 
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A Schrodinger equation may be transformed by unitary operators into dynamical equations in 
different interaction pictures which share with it a common physical frame, i.e., the same underlying 
interactions, processes and dynamics. In contrast to this standard scenario, other relations are 
also possible, such as a common interaction-picture dynamical equation corresponding to several 
Schrodinger equations that represent different physics. This may enable us to design alternative 
and feasible experimental routes for operations that are a priori difficult or impossible to perform. 
The power of this concept is exemplified by engineering Hamiltonians that improve the performance 
or make realizable several shortcuts to adiabaticity. 
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Introduction. — Schrodinger, Interaction, or Heisen- 
berg "representations" or "pictures" of a quantum sys- 
tem are linked to each other by unitary transformations 
that guarantee their formal equivalence. Changing the 
picture may be viewed as a change of basis, so in prin- 
ciple the same information can be extracted from any of 
them. This fundamental equivalence is compatible with 
distinguishing features, for example with respect to their 
usefulness {i) to calculate, approximate, or modify the 
dynamics, and (ii) to describe the system, its properties, 
and their closeness or otherwise to common language or 
classical notions. The Schrodinger picture (SP) is often 
privileged as the primary description, representative of 
the physical or experimental setting, whereas the mul- 
tiple interaction pictures (IP) have the connotation of 
auxiliary mathematical constructs to facilitate the calcu- 
lations. The standard relation among them is schemati- 
cally depicted in Fig. 1(a), where each node may repre- 
sent the dynamical equations (DE) for state vectors, the 
Hamiltonians, or the state vectors themselves. The exter- 
nal box means that they all represent the same common 
underlying physics: the same interactions and external 
forces, and the same system dynamics ("physical frame" 
for short). 

In sharp contrast to the common wisdom scenario we 
have just described, we propose in this letter alternative 
relations such as the ones in Figs. 1(b), 1(c) or more 
complex combinations, where the nodes may belong to 
different physical frames, as a way to design alternative 
and feasible experimental routes for operations that are 
a priori difficult or impossible to perform. This will be 
applied in particular to engineer Hamiltonians which im- 
prove or make feasible "shortcuts to adiabaticity" [l[. 
Adiabatic processes are very common and useful in labo- 
ratories, but their intrinsic slowness imposes limitations. 
The design of alternative fast routes is an active research 
field of interest in cold atom physics, nuclear magnetic 
resonance, quantum information processing, and beyond 
the quantum domain, e.g. to couple different devices in 
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FIG. 1: (color online). Schematic relation between differ- 
ent Schrodinger and interaction picture dynamical equations. 
Each node may also represent the Hamiltonians or the states. 
The rectangular boxes (frames) enclose nodes that represent 
the same underlying physics. The solid lines are unitary re- 
lations for the linked states, and the dashed lines represent a 
non-unitary addition of a term to the Hamiltonian. 



optical communications 0- Shortcuts designed on paper, 
may or may not be feasible in practice, so the possibility 
to generate multiple physical frames proves useful. Our 
first benchmark problem is the acceleration of adiabatic 
population inversion in a two-level system, in itself a phe- 
nomenon of broad interest from NMR applications [3| to 
quantum information Q . We shall point out several tech- 
niques to eliminate undesired Hamiltonian terms, possi- 
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ble improvements to reinterpreted experiments, and also 
feasible alternatives to inversion schemes that required 
cumbersome level-shift engineering or multiple fields 
We shall finally show how to overcome the difficulties 
to implement "counterdiabatic" terms and perform fast 
expansions or compressions of cold atoms without final 
excitation foi]. 

Multiple Schrddinger pictures. — To explain the multi- 
frame schemes of Figs. 1(b) and 1(c) we need first to 
review some basic equations and notation. SP and IP 
states are related by a unitary transformation, ji/)/) = 
U^ips), IV's) = U\iIji), and evolve according to 

ihdt\4's)=H\^Js), ihdt\^i)=Hi\i,i), (1) 

where H is the Hamiltonian in the SP and 

Hi = U\H - K)U, with K = ihUU\ (2) 

the corresponding IP Hamiltonian. Note that H, Hj and 
K may be generally time-dependent. Similar relations 
hold for a unitary operator U' which defines an interac- 
tion picture /'. 

As in Fig. 1(b), an interaction picture DE may be 
related unitarily to two or more Schrodinger equations 
with Hamiltonians that represent different experimental 
settings and external interactions. There is no contra- 
diction with the equivalence discussed above when we 
pay attention, not only to the DEs but to the observ- 
ables as well. A given picture is fully characterized by 
both the DE and the operators for the observables. Thus 
Fig. 1(b) admits several interpretations depending on the 
treatment given to the observables: If the observables are 
transformed, from Bs in S to Bs' = U'WBsUU'^ in 5', 
we will get the same expectation values from two differ- 
ent systems and dynamics performing in general rather 
different measurement operations. If instead one runs the 
same measurements in S and S" on the same, untrans- 
formed observables Bs, the expectation values will differ 
in general, and their IP representatives would be either 
WBsU or U'^'BsU', sharing in any case a common IP- 
state dynamics. It may also occur that, having applied 
the transformations, the equality Bs' = Bs holds for spe- 
cific (but not for all) observables, as in several examples 
below. 

Another multi- frame scheme is depicted in Fig. 1(c). 
In the upper box a Schrodinger node is related unitar- 
ily to a first IP-node, Ji, linked also unitarily to a sec- 
ond one, J2. The two consecutive IP Hamiltonians may 
be modified or perturbed, e.g. by the addition of some 
terms (dashed lines). This changes the physics into /' 
and /", each in a different physical frame (middle and 
lower horizontal boxes) with corresponding Schrodinger 
dynamics. In the example below the consecutive IPs are 
generated by means of adiabatic and superadiabatic it- 
erations 0, Q , and the addition of a "counterdiabatic" 
term in the Hamiltonian is performed so as to avoid tran- 
sitions, canceling out the K in Eq.® J.]. This enables 



us to accelerate slow processes without inducing any final 
excitation. 

Superadiabatic iterations and counterdiabatic 
corrections. — Our starting model Hamiltonian is 



Zj{t) Xj{t)-iYj{t) 
X,it)+2Y,{t) -Z,{t) 



(3) 



for j = and Ig = Oj i-e., Hq = X^Ux + Zqu^ in terms 
of Pauli matrices. It could represent several physical sys- 
tems such as a spin in a magnetic field, a two-level atom, 
or a condensate in the bands of an accelerated optical 
lattice In the later case, Xq may be controlled 

by the trap depth, Zq by the lattice acceleration [13] , and 
a lo component could in. principle be implemented by a 
second shifted lattice 
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The index j will be used later 
to define a series of IP Hamiltonians in successive itera- 
tions. The Hamiltonian evolution or "trajectory" is here 
specified by the Cartesian coordinates Xj, Yj, Zj. Later 
we shall also use the corresponding polar, azimuthal, and 
radial spherical coordinates, Qj, <I>j, and Rj. The radius 
may or may not be constant with respect to time, so the 
trajectory is not generally on a sphere. The Hamiltonian 
matrices are expressed in the "bare basis" of the two- level 
system, |1) = ( ^ ), \2) = {" ). We assume that this is 
also the eigenbasis of Hj{t), or very close to it for com- 
putational purposes, at initial and final times t — and 

Focusing by now on j — 0, an adiabatic population 
inversion is achieved with Hq by varying slowly Xq and 
Zq so that the resonance is crossed at Zq = 0, and the 
eigenvectors of Hq interchange their character. Differ- 
ent schemes, such as Landau-Zener (LZ), AUen-Eberly 
[l^ . and others, may be followed to specify the time- 
dependences. 

The first IP that we shall consider depends on the adi- 
abatic basis {l'T-o(O)} that diagonalizes i?o(t) and keeps 
K{t) non-diagonal. Specifically we use U = Aq, 



Ao{t) = J2 \Mt)){no{0)l 



(4) 



where we assume that the |no(0)) coincide with the bare 
basis. The corresponding K operator in Eq. (21) is de- 
noted as Kq. Also, \ipi^) = AJIV's) and Hi = AI{Ho - 
Kq)Ao. Constructing {|f^o(^))} requires a proper choice 
of phases. From an arbitrary adiabatic basis {|fia(0)} 
that diagonalizes Hq with eigenvalues En{t), 



\Mt)) 



na[t)), 



(5) 



where 7„ = i dt' (na{t')\ha{t')} is the geometric phase 
and the dot denotes time derivative. This phase choice 
is privileged since (?^o(^)|"o(^)) =0- I* also makes Kq 
non-diagonal and minimizes its norm 0, [l3t . In terms of 
the polar angle 0o, 



Kq = ihAoAl = h[Qo/2)ay. 



(6) 
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The adiabatic approximation neglects Kq in the IP 
Hamiltonian Hi to trivially solve 



(7) 



an uncoupled system in the bare basis. Alternatively one 
may add AIKqAq to Hi 0, i, The effect is to 

cancel any coupling so that Eq. (O becomes exact rather 
than an approximation. In the corresponding SP S", see 
the middle box in Fig. 1(c), this amounts to adding 

the counterdiabatic term H^^^ -.^ Kq to Ho- Hq + Kq 
preserves the populations of the approximate adiabatic 
dynamics even for short process times. 

In a new iteration, and similarly for higher orders, we 
write Hi = AI{Hq - Kq)Ao in the form of Eq. ©, 
j = 1, and diagonalize it to produce a "superadiabatic" 
basis {l'^i(i)}}: the transformation 



Ai^J2 \ni{t)){nim 



(8) 



n=l,2 



As before we assume that this basis coincides at the 
boundary times with the bare basis and that Ki = 
ifiAiAl is non diagonal in {|ni(t))}. Ai produces a new 
IP (/2 in Fig. 1(c)) with I^A/J = ^IlV^/i), and Hamilto- 
nian H2 = A\{Hi^Ki)Ai. Ki can be either neglected 
to produce a superadiabatic approximation, or canceled 
by adding a counterdiabatic term. In the corresponding 
SP [S" in Fig. 1(c)) the Hamiltonian becomes H^ + H^^J, 
where H^^^ = AqKiAI [g. In that manner a differ- 
ent shortcut Hamiltonian is created. For our reference 
Hamiltonian Hq with Yq = 0, and using polar angles for 
-ffo and Hi, 



Hcd - 7»(ei/2)(coseoa, - sineoa,). 



(9) 



if 00 < 0. The absence of a Yay component, like in 
Hq and unlike Hj,^\ is in some applications a practical 
advantage. For example, in an optical lattice implemen- 
tation of the two-level system only one optical lattice is 
required [Tlj : and in a two-level atom realization, dis- 
cussed below, it avoids the application of a second laser 
Q . One more advantage of the superadiabatic shortcut is 
that H^^ is less intense (it has a smaller norm) than H^^ 
[^. Alternative eliminations of ay are discussed next. 

Z-axis rotation. — Starting from the SP dynamical 
equation with Hq + H^^ that we write now in the form 



Hn 



rr(0) _ 
^cd - 



Zo 



-Zo 



(10) 



where (j) — arctan(?i0o/2Xo), < < 27r, and P ~ 
[Xn -I- (?i6)o/2)^]^/^, we may apply the transformation 







(11) 



which amounts to a rotation about the Z axis by 0. No- 
tice that because TJz is diagonal in the bare basis, the 
bare-state populations do not change from the SP to the 
IP. In the corresponding IP, and with = ihUzUl, the 
interaction Hamiltonian becomes 



Ul{HQ + Hf> -Kz)Uz = 



Zq - h(t>/2 

p 



p 



without Yay component. It can be realized directly 
in the laboratory and we may treat it as well as a SP 
Hamiltonian linked to the /' Hamiltonian AqHqAq, a 
common IP node for the two SP Hamiltonians in Eqs. 
([TU]) and p^ . connected via Aq and UIAq respectively. 
This Hamiltonian trio and the corresponding dynamical 
equations constitute a neat example of the dual physical 
frame scheme of Fig. 1(b). Eq. provides an alterna- 
tive shortcut path, that guarantees the same bare-state 
populations than Hq + H^fJ , and indeed it has been im- 
plemented experimentally for a condensate on an accel- 
erating lattice ll|, to avoid the realization of a ay term 
with a second optical lattice. The transition from Eq. 
([TU]) to ([T^ was justified based on properties specific to 
the optical lattice setting in [ll]. In fact the elimina- 
tion of ay in the Hamiltonian can be done formally for 
any physical realization, and its usefulness will depend 
on the feasability to implement the modified X and Z 
terms, demonstrated for a condensate on an accelerat- 
ing lattice, but more involved for a two-level atom in an 
oscillating field, see below. 

The approach based on a Z-rotation is compared in 
Figs. 2 and 3 with the one based on adding to Hq the 
counterdiabatic term i?^^'* with the LZ scheme for Hq 
(i.e. a constant Xq and a linear in time Zq). Fig. [2] 
shows the Hamiltonian matrix elements and Fig. [3] the 
populations of |1). Hq and the corresponding population 
are also shown as a reference. The process time is chosen 
to be short so that adiabaticity and population inver- 
sion fail for this Hamiltonian. Instead, the two shortcuts 
lead to perfect population inversion. Their Z compo- 
nents are similar, but the X components have a rather 
different structure. A possible advantage of the superadi- 
abatic-|-counterdiabatic approach using Hq + h'^^J is the 
smaller value of the X-maximum, which reduces ampli- 
tude noise and the field intensity. 

Two-level atoms. — In quantum optics, Eq. ([3]), with 
Xq = hVL}i/2, Yq = 0, and Zq — —fiA/2, represents a 
rotating frame IP Hamiltonian for a two-level atom in 
an oscillating field with angular frequency uj{t) = ojq — 
A{t), where ujq is the (angular) transition frequency, flji 
the (on-resonance) Rabi frequency and A the detuning, 
after having applied the electric dipole and rotating wave 
approximations (RWA). 

For K ^ Kl = -[nuj[t)/2]a^ and U = Ul ^ 
exp[— (z/?i) j^^ KL{t')dt'] (As Ul is diagonal in the bare 
basis, the populations are the same in interaction and 
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t (T/2) 



FIG. 2: (color online). Z and X Hamiltonian compo- 
nents for iifo (dashed lines); z-axis rotation Hamiltonian 
UI{Hq + H^fJ)U^ (thick blue solid lines); superadia- 
batic+counterdiabatic method Hamiltonian Ho + H'^J^ (thin 
red solid line). Zo{t) = a{t-T/2), and a = -10, T = 20/|a|, 
in units h — 1, Xo = 1. 




t (T/2) 



FIG. 3: (color online). Populations of the bare state |1) cor- 
responding to the Hamiltonians depicted in Fig. [51 Same line 
codes as in that figure. Time in units of T/2. 



Schrodinger pictures.) 
Hamiltonian is 



the corresponding S (RWA) 



(13) 



where 9{t) = j^u}{t')dt' . We can read from it the time 
dependent intensity, proportional to r2|j, the frequency of 
the field that has to be applied, 9/{2it), and the atomic 
transition frequency ujf)/{2'K). If we start instead with 
Hq + nf] in the IP DE and apply the same transforma- 
tion as before, we get the SP Hamiltonian 



Ul{Ho + 
-( 

2 \ {Hr 



+ ieo)e- 



UJQ 



,(14) 



which requires in principle two fields dephased by 7r/2 
sharing a common time-dependent frequency but with 
different time-dependent intensities Q. A rotation 
may also be applied but realizing the result is now com- 
plicated due to the time-dependence of the diagonal com- 
ponents, i.e. of the transition frequency. In this physi- 
cal context that dependence would imply time dependent 
level-shift engineering with an additional laser. It is thus 
advisable to find alternative, simpler realizations of the 
shortcuts. 

Hq + HfJ is not the only Hamiltonian that drives the 
populations along the ones of the adiabatic approxima- 
tion for . There is a whole family of them using differ- 
ent phases for the adiabatic base states, the simplest one 

'Kq)Aq = 



r(0) 



being H^^ itself Note that AI(h[^ 

so the state does not move at all in the corresponding IP, 
whereas in the DE driven by H^^J the populations will 
follow the ones for the adiabatic dynamics of iJo. By con- 
trast, alone is not enough to take the system along 
the superadiabatic path defined by Hi. We could still 
get rid of Hq and use as a shortcut to superadiabaticity 
Hfd^ - ^c? + H^cJ- To see why, use Um = AqA^, and 
notice that Uli{Hf^'^ — iT%U{)iUli)UQi = 0. In any case 



^cd^ is not so interesting for the present application as 



it combines the three Cartesian components. 

Let us now take HfJ — Kq as the reference IP Hamil- 
tonian and try to implement it with different physical 
fields as in Fig. 1(b). Applying Ul we get the SP Hamil- 
tonian 

Bs = K,^ U,K,UI = - ^^JJ, ^^^^ j , (15) 

which is indeed problematic to realize because the atomic 
transition frequency should be time dependent. In other 
words, a simple IP Hamiltonian does not necessarily 
imply a simple experiment. To remedy this, keeping 
the same simple IP DE, we may use instead U' = 
^-(i/h)!*K'(t')dt' ^ with K' = -{n/2)ujQa^. This choice 
implies now a simple resonant interaction with constant 
frequency wq and S' Hamiltonian 



Hs' = K' + U'KqU'^ 



LOq 



(16) 

Other single laser implementations may also be devel- 
oped by starting instead with Hq + H^^ . The term H^,'^ 
modifies the detuning and Rabi frequency so that the 
transformation Ul would lead to an SP Hamiltonian with 
the same structure as Eq. ([T^ . but with modified laser 
and Rabi frequencies. A further alternative to the su- 
peradiabatic iterations is the "invariants-based inverse 
engineering approach" (l8l |. 

Trap expansions. — The final example is a fast har- 
monic trap expansion, or compression, which is receiv- 
ing much attention because of fundamental and practi- 
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cal implications [1|, Il9l - l26j . The reference Hamiltonian is 
Hh = /{2m) + where w = Cj{t) is the time 

dependent angular frequency, m the particle mass, and q 
and p are position and momentum operators. The cor- 
responding counterdiabatic term to avoid excitations is 
^cd ~ ~(m + 9P)'^/(4w), whose direct laboratory 
implementation is problematic and was left as an open 
question 0]. This difficulty is overcome by the trans- 
formation Uq = exp(i^^g^), which eliminates the cross 
terms; it produces from Hs = Hh + HfJ the IP Hamilto- 
nian Hi = Ul{Hs - inUqUl)Uq = p2/(2m) + 
where 



3£^ 

4^ 



2Cj 



1/2 



(17) 



This Hamiltonian can actually be realized directly [20|, 
2l| and considered in a different physical frame as an 
ordinary harmonic oscillator with modified frequency. To 
satisfy the scheme of Fig. 1(b) wc may apply [/ = 1 and 
regard it as a Schrodinger Hamiltonian S", namely Hi — 
Hs' ■ It indeed provides a shortcut with the following 
properties: starting with a common state at time t = 0, 
the spatial densities driven by Hh + HfJ and Hs' are 
identical. In fact, by imposing a)(t/) = w(i/) = the 
final state is also equal for both dynamics, even in phase, 
and the final vibrational state populations coincide with 
those of a slow adiabatic process. 

Discussion. — We have first proposed schemes for 
which different interaction and Schrodinger picture dy- 
namical equations represent different physical processes 
and interactions. These schemes have been later com- 
bined and exemplified to produce better, realizable short- 
cuts to adiabaticity for population inversion protocols, 
and for expansions and compressions. Similar manipu- 
lations may be applied as well to facilitate or improve 
shortcuts to adiabaticity for other operations such as 
controlled atomic transport 27-2^. In fact the idea of 
designing the pictures to generate alternative, easier to 
handle physics, is applicable to a plethora of quantum 
systems, in particular, in the realms of quantum simula- 
tions, quantum control, or quantum information, where 
developing techniques to externally drive the systems for 
specific goals is a central objective. 
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